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Abstract 

_ The high-temperature expansion coefficients of the ordinary and the higher susceptibilities of 

the spin- 1/2 nearest-neighbor Ising model are calculated exactly up to the 20th order for a general 
d-dimensional (hyper)-simple-cubical lattice. These series are analyzed to study the dependence 
of critical parameters on the lattice dimensionality. Using the general d expression of the ordinary 
susceptibility, we have more than doubled the length of the existing series expansion of the critical 
temperature in powers of 1/d. 
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I. INTRODUCTION 



We have derived high-temperature(HT) expansions of the ordinary and the higher sus- 
ceptibilities (see the definitions in Section II) of the spin-1/2 Ising model exactly up to the 
20th order for the general (i-dimensional (hyper)-simple-cubical(hsc) lattices. 

These expressions are obtained by interpolation of the HT series expansion coefficients of 
the susceptibilities over the integer values of the lattice dimensionality, and not by analytic 
continuation. Thus they have no obvious uniqueness properties when d is allowed to take 
noninteger values. A priori a different dependence of physical quantities on the dimensional- 
ity might result from different possible interpolations, such as that obtained by formulating 
the Ising model on a fractal lattice^, whose Hausdorff dimension can be varied continuously, 
or also in other ways 3 . An analogous situation is known to occur for the n- vector model, 
whose HT series coefficients can be interpolated 4 by rational functions of n. Of course, 
non- integer values of d (or similarly of n) might in some cases have no physical meaning 5 . 

The expansions^"- of the physical quantities in powers of 1/d, i.e. around the mean-field 
(MF) approximation (or in powers^ of 1/n, i.e. the expansions around the spherical model 
limit), are related with these analytic representations in terms of d (or of n). 

Our results not only provide reference data in a compact form for the higher susceptibil- 
ities, which are generally difficult to compute by methods different from series expansions, 
but also make a variety of other investigations possible. For example, in discussing^"— 
how the finite-size-scaling behavior— changes for systems above the upper critical dimen- 
sionality, accurate estimates of the critical parameters are needed as benchmarks. Our 
data can also help to assess the accuracy of estimates of physical parameters obtained from 
approximations of a different nature, such as the the e = 4 — d expansion 9,18 of the renormal- 
ization group theory, the fixed- dimension renormalization group 1 ^, the 1/d-expansion^, the 
MonteCarlo(MC) simulations, etc.. It is appropriate at this point to observe that the MC 
simulations become increasingly time- and memory- demanding as the lattice dimensionality 
d grows, whereas the non-analytic corrections to scaling in the asymptotic critical behavior 
of physical quantities, that usually make the HT series analyses a delicate matter and are 
the main source of their uncertainties, become simpler and smaller— with increasing d. Thus 
even moderately long HT expansions can lead to very accurate estimates already for not 
very large d. 

It is also worth mentioning that extremely long, although approximate HT expansions 
of the ordinary susceptibility have been recently generated^ 1 - by a MC method for the hsc 
lattices of dimensionalities d — 5, .., 8 and used to test the accuracy of the "extended scaling" 
ideas22r— above the upper critical dimension. The results of this investigation can now be 
compared with those from the analysis of our far shorter, but exact expansions. 

The expressions presented here are obtained from recently derived^ 5 ."— HT and low-field 
series expansions of the magnetization in presence of an external magnetic field for the spin- 
1/2 Ising model with nearest-neighbor interactions. Actually, we have produced a wider 
ranging set of data including as well other spin-systems in the Ising model universality class, 
such as the general spin-s Ising model and the lattice scalar-field theories with polynomial 
self-interaction and thus, also for these models we are able to write exact expressions valid 
for general d- dimensional hsc lattices. 

Our derivation of the HT and low-field series, that significantly extend the longest known 
results, even in zero field 4 - 1 ^"— , have been performed for the lattice dimensionalities d = 
1,2, ...10. The expansions are carried to the 24th order in the case of the (hyper)-body- 



2 



centered-cubical lattices, but they are slightly shorter in the case of the hsc lattices: we 
have obtained the 24th order for d < 5, the 22nd for d = 5, the 21st for d = 6, the 20th for 
7 < d < 10. 

The layout of the paper is the following: in the second Section, we set our notations and 
tabulate the expansion coefficients of the ordinary susceptibility as closed-form polynomials 
in the coordination number q = 2d. The corresponding data for the higher susceptibilities 
are reported in the appendix. The printed version of this paper— contains no appendix and 
the data for the higher susceptibilities are retrievable as Supplemental Material at the URL: 
http://link.aps.org/supplemental/10.1103/PhysRevE.86.011139, In the third Section, we 
discuss the results of some analyses of these series. The last Section contains our conclusions. 



II. THE ISING MODEL IN GENERAL DIMENSION 

The partition function of a spin-1/2 Ising system with nearest-neighbor interactions, in 
the presence of an external magnetic field H, on a finite (i-dimensional lattice of N sites can 
be written as 

Z N (H,T; d) = ^exp J/k B T ^ SjSj + H/ksT^sA . (1) 

conf <ij> i 

Here Sj = ±1 denotes an Ising spin variable associated to the lattice site i. The first 
sum extends to all configurations of the spins, the second to all distinct pairs < ij > of 
nearest-neighbor spins and the third to all spins. We shall set K = J/ksT, with T the 
temperature, J the exchange coupling, k B the Boltzmann constant, and h = H/ksT the 
reduced magnetic field. In terms of the variable v = tannic, the HT expansion coefficients 
are simple integers and so this variable is more convenient for the data tabulation. 

In the thermodynamic limit N — > oo, the specific free-energy J-(h, K; d) is defined by 

_ —f(h,K;d) = lim —\nZ N (h, K; d) (2) 

J N^oo N 

The specific magnetization A4(h,K; d) is defined by 

a^/l j\ K dT{h,K]d) . . 

M(h,K;d) = -j y — >- (3) 

Our HT calculation of the field-dependent magnetization, has yielded significant extensions 
of the existing HT expansions in zero field for the 2p-spin connected-correlation-functions at 
zero wave-number and zero field X2p{K;d) (usually called "higher susceptibilities"). These 
quantities are defined by the successive field-derivatives of the specific magnetization 

X2p (K; d) = (d 2 ^ l M(h, K- d)/dh 2 v~ l ) h=0 = < s i s 2-S2 P > c • (4) 

S2,S3,...,S2p 

at zero field. 

The even field-derivatives of the magnetization vanish at zero field in the symmetric HT 
phase, while all derivatives are nontrivial in the broken-symmetry low-temperature phase. 

The HT expansion coefficients of the susceptibilities at a given order I in K are usually 
computed as sums of contributions classified in terms of graphs having I edges. To each graph 
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we associate a weight depending on the symmetry of the graph and on its free-multiplicity, 
i.e. the number of distinct ways (per lattice site) in which the graph can be embedded in 
the lattice, associating each vertex to a site and each line to a nearest-neighbor bond 3 ^. 
Only the latter quantity, technically denoted as "free-lattice- embedding number" within the 
linked-cluster HT expansion method, depends on the lattice dimensionality. An analysis of 
these numbers for the various classes of graphs, like that performed in Refs. leads to the 
conclusion that, at any given expansion order /, the HT series coefficients of the ordinary and 
the higher susceptibilities can be written as simple polynomials in the lattice dimensionality 
d of degree at most I. It is then clear that a straightforward prescription to represent 
the HT series coefficients of a susceptibility as polynomials in d up to the order K lmax , 
consists in repeating the computation of the HT series for lattices of dimensionalities d = 
0, 1, l m ax and then in interpolating each series coefficient with respect to d. Unfortunately, 
this straightforward strategy works only for relatively small orders of expansion, since, in 
the case of the hsc lattice the combinatorial complexity of the computation of the graph- 
embedding numbers grows large exponentially with the dimensionality. 

Here we take advantage of a well known result that is helpful in mitigating this difficulty. 
It was shown long age 35 that the HT expansions of the successive derivatives of the magnetic 
field with respect to the magnetization d 2p+l h/ dM. 2p+1 at zero magnetization, for p = 0, 1... 
contain only star graphs, i.e. connected graphs having no articulation vertex. This property 
can be very simply understood observing that in field-theoretic language these derivatives 
define the connected amputated one-particle-irreducible correlations at zero wave-number. 

This property was sometimes used to restrict the number of graphs contributing to the 
HT expansion of higher susceptibilities. What is interesting for our aims is the fact that the 
lattice-embedding number of a star graph with I edges is a polynomial in d of degree [1/2] 
at most. Here [1/2] denotes the integer part of 1/2. The higher susceptibilities are simply 
related to the quantities d 2p+1 h/dM 2p+1 : 



8M y ' ' X 2(K;d) 

'V;^) = ^#4 (6) 



&M 3V ' 1 X 2(K;d) 

d 5 h X 6(K; d) x^K^df 

-\K;d)= ttt-^b -IO - TTTTlvr (') 



dM^ ' y X2(K;dy X2(K;d)< 
97k m d) = - 56 ^(K;d)MK;d) + 2m MK;df (g) 



dM 7X ' ' X2{K-df X2{K;df X2(K;dy° 



d 9 h _ XiojK; d) X4(if; d)xs(K; d) xe{K; d) 2 

dM* [K ' d) X 2(K;dyo X 2(K;d)^ ^(Kjd) 11 

■ A ^ n X4(K-d) 2 X e(K;d) „ nn X^df 
+462 ° X2(K;dY 2 1540 °^P (9) 



dUh -(K- d) = Xl2 ^ d ^ - 22Q X4(K;d)x w (K;d) _ ^Xe^ d) X s(K; d) 



dM^ y ' ' X2{K-d) 12 X2(K;dyz - X 2(K-d) 
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+17160 



X4 (K; d) 2 X s(K;d) 



+ 36036 



X2(K;dy± 



- 560560 



X^K-dfx^K-d) 
X2{K-d)^ 




(10) 



and so on. 



Then it transpires that it is sufficient to interpolate the HT series expansion coefficients 
of the combinations of (higher) susceptibilities on the rhs of eqs.(15])..... (fTD|) etc. only over 
the dimensionalities 1 < d < [1/2], to obtain the representations of these quantities for 
general d through the l-th order in K. Finally, from these results the representations of the 
single higher susceptibilities can be easily recovered. This simple remark leads to a decisive 
reduction of the combinatorial complexity of the necessary calculations. In our case, only 
the knowledge of the HT expansions of the (higher) susceptibilities for all hsc lattices with 
d < 10 up to the 20th order is sufficient to obtain the expression of these susceptibilities for 
general d up to the same order. The fact that the coefficient of order K l in X2p{K; d) must 
be a polynomial in d of degree /, while the corresponding coefficient at the same order in the 
expansion of d 2p ~ 1 h/dA4 2p ~ 1 is a polynomial of degree [1/2], provides a simple consistency 
check of our computations. 

A brief technical comment on a detail of the calculation is appropriate at this point, since 
most of the computing time goes into counting the number of lattice embeddings of each 
graph for relatively large lattice dimensionality. The first step of the computation consists 
in ordering appropriately the vertices of the graph and in placing the first of them at the 
lattice origin. The second step consists in counting the possible positions, of coordinates 
(xx, ..,Xd), for the second vertex. It is crucial to optimize this step by using the hypercubical 
symmetry to restrict the possible positions of the second vertex to the fundamental region 
%i > x 2 > > Xd > 0. Separating each > case into a > and an =, one gets 2 d cases. 
The case with all > corresponds to the inside of the fundamental region, whose sites are 
representatives of a group orbit of length d\2 d . A small program precomputes the length of 
the orbits for the each of the 2 d cases, which is then used in the embedding program. After 
fixing the first two points of the embedding, the possible positions of the remaining vertices 
are restricted to relatively few sites by the constraints given by the distances from the first 
two points. 

The timings for computing the HT expansion of the (^-dimensional Ising model at order 
I have been roughly 0(5.5 z 2.5 d ). In particular, the 10-dimensional Ising model at order 20 
took 42 days of single-core time on a quad-core computer with a CPU-clock frequency of 
2.8GHz. 

III. THE HT EXPANSIONS 

Expressions like those obtained here have already appeared at lower orders in the litera- 
ture for the spin- 1/2 Ising model^ 1 ^ - — with nearest-neighbor interaction. In the case of the 
ordinary susceptibility Xi{K\ d), they reached at most^ the 15th order, while for Xi(K; d) 
and Xe{K; d) they were computed^ only up to the 11th order. 

For brevity, we shall report our results in the Table H] only for X2{K;d) in the case 
of the spin-1/2 Ising model. The corresponding data for the higher susceptibilities 
Xi{K; d),...Xi&{K; d) are reported in the Appendix. 
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We shall tabulate here also the coefficients of order < 16, although they reproduce those 
already listed in Ref. 37]]. A simplification of the expressions of the coefficients is obtained 
by using the variable q = 2d rather than the variable d used in Ref. [37[ • As already pointed 
out, we can produce analogous formulas also for the other models in the Ising universality 
class for which we have computed the HT expansions of the magnetization, but they will 
not be presented here. 



(2) 

TABLE I: The coefficients <£ J (d) of the HT series expan- 
sion in powers of v = tanh-K" for the ordinary susceptibility 
X2(v]d) = ^2^Lq Cn (d)v n of the spin- 1/2 Ising model with 
nearest-neighbor interaction, on a (hyper)-simple-cubical lat- 
tice of general dimensionality d. For convenience, we have re- 

(2) 

produced here also the expressions of the coefficients dk (d) 
with n < 16 that were already tabulated in Ref. [371] in terms 
of the variable d, whereas here we have used the variable 
q = 2d. It should also be stressed that also in Ref. |37l] the 



series coefficients refer to the expansion variable v like in this 
Table and not to the variable K as erroneously stated. 



4 2) = q 3 - 2q 2 + q 

C W = g 4 - 3q 3 + 2q 2 + q 

4 2) = q 5 - 4g 4 + Aq 3 + 2q 2 - 3q 

4 2) = q 6 - 5q 5 + 7<? 4 - 2q 3 + I2q 2 - 23q 

4 2) = q 7 - 6<? 6 + Uq 5 - 8q 4 + 33g 3 - 88q 2 + 61q 

c {2) =q s - 7q 7 + 16q 6 - 17q 5 + 33<? 4 + 77q 3 - 652q 2 + 813q 

c {2) =q 9 - 8q 8 + 22q 7 - 30q 6 + 41q 5 + 782/3<? 4 - 5576/3g 3 + 10594/3g 2 - 6145/3g 
cg> = gio _ 9^9 + 29 g 8 _ 48q 7 + 6Qq e + 595/3^ + 1132/3^ _ 4i404/3g 3 + 141125/3<? 2 
-45837g 

cj 2) = q 11 - 10q w + 37q 9 - 72q 8 + 94q 7 + 128g 6 + 2791g 5 - 112408/3g 4 + 140651g 3 
-628610/3^ 2 + 105297g 

eg) = q 12 - Uq 11 + 46q w - 103g 9 + 148<? 8 + 110/3g 7 + 7411/3g 6 - 4413g 5 
-788360/3<? 4 + 5431963/3<? 3 - 13319981/3g 2 + 3721963g 

C J 2) = q ls - I2q 12 + 56<? n - 142g 10 + 228g 9 - 278/3<? 8 + 6760/3g 7 + 90088/3<? 6 

-2268283/3g 5 + 4828108<? 4 - 13659106<? 3 + 52372552/3g 2 - 7921783g 

C J 2) = q u - 13q 13 + 67<? 12 - 190g n + 341g 10 - 282g 9 + 2163g 8 + 80084/3<? 7 - 633892/3<? 6 

-14257841/3g 5 + 179124608/3g 4 - 269431329g 3 + 1643435774/3<? 2 - 414746143<? 

C J 2) = q 15 - Uq u + 79q 13 - 248q 12 + 495g 41 - 1678/3<? 10 + 6724/3g 9 + 23674g 8 

+347769<? 7 - 78796706/5<? 6 + 468273229/3<? 5 - 2167816142/3<? 4 + 5176873061/3<? 3 

-9973026754/5<? 2 +2545796651/3g 

cj 2) (d) = q 16 - 15q 15 + 92q 14 - 317g 13 + 699q 12 - 2879/3q n + 7663/3g 10 

Continued on next page 
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(2) 

TABLE I: The coefficients c£ J (d) of the HT series expan- 
sion in powers of v = tanhi<C for the ordinary susceptibility 
X2(v;d) = J2™=o Cn \d)v n of the spin- 1/2 Ising model with 
nearest-neighbor interaction, on a (hyper)-simple-cubical lat- 
tice of general dimensionality d. For convenience, we have re- 

(2) 

produced here also the expressions of the coefficients c n (d) 



with n < 16 that were already tabulated in Ref. 37J in terms 
of the variable d, whereas here we have used the variable 
q = 2d. It should also be stressed that also in Ref.ji3] the 
series coefficients refer to the expansion variable v like in this 
Table and not to the variable K as erroneously stated. 

+62404/3g 9 + 951902/3<? 8 - 29340047/5<? 7 - 1222190629/15<? 6 + 1842802906<? 5 
-40303287247/3g 4 + 727440333881/15^ 3 - 436050363522/5<? 2 + 61362596609.? 
c V)(d) = q 17 - 16q ie + 106<? 15 - 398g 14 + 963<? 13 - 1526<? 12 + 3198g n + 53594/3<? 10 
+872213/3<? 9 + 62822606/15<? 8 - 1721035544/5<? 7 + 24837980998/5g 6 
-103041998423/3<? 5 +658765700268/5^ - 1406657809016/5<? 3 
+1518632487582/5<? 2 - 124150140027^ 

c^(d) = q 18 - 17q 17 + 121g 16 - 492g 15 + 1298<? 14 - 6931/3g 13 + 12916/3g 12 
+14580g n + 803611/3<? 10 + 58164043/15g 9 - 6630346477/45g 8 - 18614162023/15<? 7 
+2490068122951/45g 6 - 9120432843283/15<? 5 + 153498549866917/45g 4 
-160125507535387/15<? 3 + 791246681603369/45g 2 -35100413743831/3^ 
c g>(d) = g i9 _ l8g 18 + 137<? 17 - 600g 16 + 1716<? 15 - 10124/3<? 14 + 18163/3<? 13 
+31678/3<? 12 + 746245/3g n + 10791224/3<? 10 + 155233339/3.? 9 
-119093401726/15g 8 + 797131628934/5<? 7 - 7738214717002/5<? 6 
+26108335469563/3g 5 - 441756642770324/15g 4 + 870753648372433/ 15<? 3 
-893908363283714/15g 2 + 23668071765201^ 

c g(d) = q 20 - 19q 19 + 154g 18 - 723g 17 + 2230g 16 - 4792g 15 + 8677g 14 
+15815/3<? 13 + 701276/3<? 12 + 50129237/15g n + 48537033<j 10 - 18335572678/5g 9 
-109608216238/9<? 8 + 8227897945731/5<? 7 - 1179114186108353/45<? 6 
+1059691462407483/5<? 5 - 9089741189180104/9.? 4 + 2855597663272273<? 3 
-65971006971414364/15<? 2 + 2798133827599029^ 



IV. HIGH DENSITY EXPANSION OF THE CRITICAL TEMPERATURE 

By solving recursively the equation 

1/X2(v c ;d) = (11) 

with the Ansatz v c (d) = tannine?) = a^jq + G^/g 2 + ... an expansion of the critical temper- 
ature in inverse powers of q was computed in Ref. p] up to the fifth order. We confirm the 
results of Ref. jsj , (which are expressed in terms of the variable K rather than the variable 
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v used above) and, taking advantage of our data, we are able to carry the 1/q expansion to 
the 12th order: 

1 1 4 13 979 2009 176749 6648736 765907148 

qK c (d) ~ ~ q ~ 3f ~ 3g3 ~ 45q^ ~ I5q^ ~ 189g 6 945<f 14175g 8 
5446232381 829271458256 164976684314 6495334834824112 
14175g 9 467775g 10 + 22275g n + 638512875g 12 "' ^ ' 

A fifth order expansion of the same kind for K c (d) was also obtained^ for the n- vector model. 
All these expansions are presumably of asymptotic character, but so far this property has 
been established^ only in the case of the spherical model, i.e. in the large n limit. 



V. SERIES ANALYSES 

We shall now very briefly discuss the numerical estimates of some non-universal critical 
parameters of the spin-1/2 Ising models for d > 4. In particular, we shall use the expansion 
of x 2 (K\ d) to locate the critical points K c (d). We shall also estimate the critical amplitudes 
of the five lowest-order susceptibilities and a few universal ratios of these. In our analysis, 
we have simply assumed that all (higher) susceptibilities show MF exponents of divergence, 
as also our recent work— ^ has contributed to confirm numerically. 

The critical parameters are defined by the asymptotic critical behaviors of the suscepti- 
bilities 

X 2 P (K; d) » A 2p (d)r(d)~^ [l + a 2p (d)r 9 ^..^ (13) 

as K -> K c (d)~. Here r(d) = (1 - K/K c (d)), 7 2p = 72 + 3(p - 1) is the MF exponent, 
(assumed to depend on the order 2p of the susceptibility, but not on the lattice dimensionality 
for d > 4) and 72 = 1. A 2p (d) denotes the amplitude of the leading singularity, a 2p (d) 
the amplitude of the leading correction to scaling and 6(d) is the exponent of the leading 
correction-to-scaling. The value of 6(d) is expected^ to be 1/2 for d — 5, while for d — 6 it 
should be 1 with a possible logarithmic multiplicative correction. Generally, for d > 6 it is 
expected that 6(d) = (d - 4)/2. 

We can only briefly outline the now standard numerical approximation techniques that we 
have used for these analyses, since a more detailed discussion was already given in Refs. (2HIJ26L 



281 ] . We have mainly employed the differential approximant(DA) method, that generalizes 39 
the elementary well known Pade approximant method, to resum the HT expansions up 
to the border of their convergence disks. This technique estimates the values of the finite 
quantities or the singularity parameters for the divergent quantities from the solution, called 
differential approximant, of an initial value problem for an ordinary linear inhomogeneous 
differential equation of the first- or of a higher-order. Various differential equations can be 
formed from a given series expansion. For each of them, the coefficients are polynomials in 
the expansion variable, defined in such a way that the series expansion of the solution of the 
equation equals, up to some appropriate order, the series to be approximated. 

Sometimes, to determine the location of the critical points, it is more convenient to use a 
smoother and faster converging modification2&22^, called modified-ratio-approximant(MRA) 
of the traditional methods of extrapolation of the series coefficient-ratio-sequence. The 
MRAs produce sequences (K^\d)) of approximations of the critical point that can be easily 
extrapolated to large orders r of expansion and therefore in some cases they may yield 



more accurate estimates than the DAs for which the analogous extrapolation is somewhat 
arbitrary. Let us finally stress that, when using the DAs the evaluation of the uncertainties 
has not the same meaning as for MCs, but remains subjective to some extent and only 
indicates a small multiple of the spread of the values of a conveniently large sample of the 
highest-order approximants, formed from all or most expansion coefficients. If the sample 
averages remain stable as the order of the series grows and it can be assumed that stability 
indicates convergence, than the spread can be seen as a reasonable measure of the uncertainty 
of the results. 

For the critical inverse temperatures K c (d) of the systems under study, we consider our 
best estimates those reported in Table [XTJ They are obtained from the HT expansion of the 
ordinary susceptibility X2(K;d) by extrapolating to large order of expansion, a few (from 

four to six) highest-order terms of the MRA sequence of estimates (Kc (d)) of the critical 
inverse temperature, basing on a fit to their simple asymptotic behavior^ 

Wl-KM- ™^ 1 -™* ^). (14) 

A small multiple of the fit error is taken as the uncertainty of the final estimate. 

In the case of six-dimensional lattices, we expect 9 — 1. Therefore the second term on the 
right-hand side of eq. (H4j) vanishes and it must be replaced by a higher-order term reflecting 
the exponent of the next-to- leading correction to scaling in eq.f lT3|) . In the d — 5 case, in 
which one expects 9 = 1/2, the coefficient of \/r l+e in eq. (|14p appears to be numerically 
negligible, so that the situation is similar to that of the six-dimensional case. In general, to 
avoid making assumptions on the values of the exponents of the next-to-leading correction 
to scaling, we have assumed an asymptotic form Kc — K c + w/r l+e and determined K c , 
w, and the "effective" exponent e = e(d) by a best fit to the few highest-order terms of the 
sequence (K^). We thus obtain the values e = 1.1(2) for d = 5 and e = 1.5(2) for d = 6. For 
d > 6, the values of e thus obtained are larger. Therefore our estimates consistently confirm 
the above mentioned expectations about 9(d), and suggest that in d — 5 and d = 6 the 
asymptotic behavior of eq. (fT4"j) is actually determined by the next-to- leading, rather than 
the leading correction to scaling. On the other hand, for d — 5 and d — 6, a measure of the 
exponent 9(d) of the leading corrections to scaling, whose amplitudes <i2 P (d) are probably not 
negligible in spite of the fact that they are not seen by the MRAs, can be performed studying 
by DAs the critical behavior of appropriate universal ratios of higher susceptibilities, such 
as those introduced below in eqs. (|16p . (|17p and (|18|) . In these ratios the dominant critical 
singularities cancel, while the leading corrections to sca ling generally survive and thus can 
be detected by DAs. This was already discussed in Ref. [261] . In conclusion, all these results 
are in reasonable agreement with the expect at ions^^i indicated above. 

In the Table [Til we have also reported a few of the most recent estimates of the critical 
inverse temperatures obtained in the literature either from shorter HT series or from MC 
simulations, for the values of d considered in our study. Since, for d > 4, no logarithmic 
factors are expected to modify the leading MF behavior of the physical quantities, our series 
analyses are likely to yield estimates of a high accuracy, that moreover seem to improve with 
increasing lattice dimensionality, both because of the decreasing influence of the corrections 
to scaling and of the increasing lattice coordination number. All the results obtained from 
MRAs are consistent, within their uncertainties, also with the analyses employing DAs. In 
d — 5 dimensions, our estimate is slightly larger than other estimates^^ of similar nominal 
accuracy, but can be considered essentially compatible with those of Refs. (2l|j37l]4lll42l | . It is 
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TABLE II: Our estimates of the critical inverse temperatures K c (d), obtained from the ordinary 
susceptibility expansions, for several hsc lattices of dimensionality d > 4. We have marked by 
an asterisk the estimates in the cases in which expansions extend beyond the 20th order—. In 
particular for d = 5 our series extend to the 22nd order, and for d = 6 to the 21st order. 



Source 


Kc(5) 


K c (6) 


Kc(7) 


K C (S) 


K c (9) 


K c (10) 


HT This work 


0.113920(1)* 


0.092298(1)* 


0.0777094(2) 


0.067155(1) 


0.059148(1) 


0.052858(1) 


HT [37.411 


0.113935(15) 


0.092295(3) 


0.077706(2) 








HT [42] 


0.113915(3) 












MC [21] 


0.113925(12) 


0.092290(5) 


0.077706(2) 


0.067144(4) 






MC [12,13] 


0.11391(?) 


0.09229(4) 










MC [14-161 




0.09229(4) 


0.0777(1) 


0.06712(4) 






MC [10] 


0.1139152(4) 












MC [111 


0.1139139(5) 













of interest to quote here also our estimate K c (5) = 0.113919(2) obtained from second-order 
quasi-diagonal DAs that use all series coefficients up to order 20 < / < 22. The same value, 
with a slightly larger uncertainty, is obtained from DAs using all series coefficients up to 
order 18 < / < 20. Generally, for higher values of d, our estimates of the critical inverse 
temperatures do not differ much from the old values, but they show a greater accuracy. For 
d > 7, the estimates of K c (d) obtained from the 1/d expansion of eq.f lT2|) reproduce our 
MRA values within the errors. 

Let us now turn to the critical amplitudes A 2p (d) of the susceptibilities X2 P (K; d) with 
n = 1,2, ...5, that can be determined, in terms of the previously estimated values of K c (d), 
by extrapolating the effective amplitudes 

A% f (K ; d) = (1 - K/K c (d)) 72p X2 P (K; d) (15) 

to K = K c (d), namely A 2p {d) = A e 2 J (K c ;d). Our analysis uses first- and second-order 
DAs of the HT expansion of A e 2 J [K; d). 



TABLE III: Our estimates of the critical amplitudes A 2p (d) of the susceptibilities X2p(K;d), 
normalized to their values AM, in the MF approximation for several hyper-simple-cubical lattices 
of dimensionality d > 4. We have marked by an asterisk the estimates obtained from series 
extending beyond the 20th order. (For d = 5 our series extend to the 22nd order, and for d = 6 
the 21st order.) 



Amplitude 


Source 


d = 5 


d = 6 


d = 7 


d = 8 


d = 9 


d = 10 


A 2 (d)/A* IF 

Md)/A™ F 
Md)/A%> F 
A 10 (d)/Af F 
Md)/A™ F 
A 2 {d)Af F 


This work 
This work 
This work 
This work 
This work 
MC [21] 
HT [201 


1.32(1)* 
1.40(1)* 
1.49(1)* 
1.57(1)* 
1.65(2)* 
1.291(3) 
1.311(9) 


1.179(2)* 
1.20(1)* 
1.22(1)* 
1.25(1)* 
1.28(1)* 

1.1606(17) 
1.168(8) 


1.124(2) 
1.138(2) 
1.147(3) 
1.165(2) 
1.19(1) 
1.1008(5) 


1.096(2) 
1.105(2) 
1.114(2) 
1.122(2) 
1.13(1) 
1.0836(5) 


1.079(1) 
1.085(1) 
1.091(2) 
1.097(3) 
1.103(4) 


1.067(3) 
1.068(3) 
1.077(3) 
1.081(4) 
1.085(4) 



In the Table IHIl we have reported our series estimates of A 2 (d), At(cf), A s (d), A$(d) 
and A w (d) normalized to their values in the MF approximation™: A^f F = 1, A± = —2, 
Af F = 40, Af F = -2240 and Af 1 / = 246400. As expected, these ratios tend to unityas 
d — > oo. For comparison, we have also reported the corresponding MC results of Ref. 21 



The results of this simulation are slightly, but systematically smaller than our HT estimates. 
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This minor disagreement cannot be due to the very small difference in the estimates of K c (d) 
used in the two calculations, but must probably be ascribed to an underestimate of the 
uncertainties inherent in the MC generation of the HT series. On the other hand, the old 
estimates^ of the same amplitudes obtained from an analysis of the 11th order HT series 
of Ref.0| are compatible with ours. Estimates for A 2 (d) have also been obtained 8 from a 
third order expansion in 1/d, but they are not accurate enough. No comparison at all was 
possible for A^(d),...,Ai (d) for which we do not know of other estimates in the literature. 

Finally, we have computed also for d > 6, the critical values of a few universal ratios 
of higher susceptibilities such as the lowest order terms in the sequences I 2r+4 (d), A 2r+4: {d) 
and B 2r+8 (d), defined by 



-2r+4 



(d) 



B 2r+8 {d) 



lim 
lim 
lim 



X2(K;d) r X2r + 4(K;d) _ A 2 (d) r A 2r+A (d) 
X 4(K;d) 2r+1 A 4 (d) 2r+1 

X2r(K; d)x2r+i{K\ d) _ A 2r (d)A 2r+4 (d) 

A 2r+2 (d) 2 

A 2r (d)A 2r+8 (d) 



X2r+2(K;d) 2 

X2r(K; d)x2r+s(K] d) 
X2r+4(K; d) 2 



A 



2r+4 



(dy 



(16) 
(17) 
(18) 



for r > 0. These universal ratios were introduced in Ref . [44J and were studied in detail for 
d = 4, 5, 6 in Ref. 021. For the first few values of r = 1,2,3, we have checked that as 
expected, also for d > 6, they take the MF values: 1+ MF = 10, 1+ MF = 280, I+ MF = 15400, 
A$ AIF = 14/5, Aiq F = 55/28 and Bf^ IF = 154 within a relative accuracy generally higher 
than 10 -3 , although the single amplitudes entering into the ratios reach their MF value only 
in the d — > oo limit. 



VI. CONCLUSIONS 

We have represented in a compact and exact form, as simple polynomials in the lattice 
dimensionality, the HT expansion coefficients of the (higher) susceptibilities in the case of 
the spin- 1/2 Ising model on the hsc lattices of general dimensionality d. Our calculations 
add five more orders to the existing expansions of the ordinary susceptibility for general 
d and nine more orders to those of the fourth- and sixth-order susceptibilities. For the 
susceptibilities of order greater than the sixth no such data already exist in the literature. 

An analysis of the series for lattice dimensionality d > 4 yields estimates of non-universal 
parameters that compare well with the previous results whenever available, but are generally 
more accurate. Our estimates of a few universal ratio amplitudes, provide a high accuracy 
check that, unsurprisingly they take MF values for d > 4. 

Finally, using the general d expression of the ordinary susceptibility, we have been able to 
expand up to the 12th order the critical temperature in powers of 1/d more than doubling 
the length of the expansion already known in the literature. 



VII. ACKNOWLEDGEMENTS 



We thank the Physics Departments of Milano-Bicocca University and of Milano University 
for hospitality and support. Partial support by the MIUR is also acknowledged. 



11 



VIII. APPENDIX 



A. The HT expansion of Xi{v;d) 



TABLE IV: The coefficients c ( n\d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility Xa{ v \ d) = 
En=0 C n\d)v n Of the spin- 1/2 Ising model on a ( hyper )- 
simple-cubical lattice of dimensionality d as polynomials in 
the coordination number q = 2d. 

cf = -8, 

4 4 ) = -2O9 2 + Uq 

cp = -40q 3 + 64q 2 - 2Aq 

c| 4) = -70g 4 + 180g 3 - 122g 2 - Uq 

cl 4) = -112<7 5 + 400q 4 - 408<? 3 - 48q 2 + 168<? 

cp = -168g 6 + 770q 5 - 1080g 4 + 152<? 3 + U6q 2 + 578<? 

A 4) = -2A0q 7 + 1344g 6 - 2440g 5 + 1120g 4 - 1040<? 3 + 5360<? 2 - 4424g 

A 4) = -330g 8 + 2184g 7 - 4914g 6 + 3850g 5 - 3888g 4 + 8464g 3 + 14374g 2 - 33774g 

J® = -440g 9 + 3360g 8 - 9072g 7 + 10016g 6 - 10144g 5 + 2864/3g 4 + 408256/3g 3 - 980528/3g 2 
+618536/3g 

c j 4) = -572g 10 + 4950g 9 - 15648g 8 + 22232g 7 - 23722g 6 - 34580/3g 5 + 763630/3g 4 - 37060/3g 3 
-6036988/3<? 2 + 2586414^ 

c [ 4) = -728<? n + 7040g 10 - 25560<? 9 + 44352g 8 - 52024g 7 - 19616g 6 + 221136g 5 + 7971968/3g 4 
-15033816g 3 + 76181488/3g 2 - 13485336g 

c [ 4) = -910g 12 + 9724g n - 39930g 10 + 81810g 9 - 107730g 8 - 18080/3g 7 + 182096/3g 6 

+17071112/3g 5 - 42028180/3g 4 - 62677128g 3 + 806882456/3<? 2 - 262544882g 

c [ 4) = -1120<? 13 + 13104<? 12 - 60104Q 11 + 142000<? 10 - 211128g 9 + 182912/3<? 8 - 668440/3<? 7 

+19780304/3g 6 + 41895984g 5 - 1655518208/3g 4 + 5642509000/3g 3 - 7772420336/3g 2 
+1218754376^ 

cf} = -1360g 14 + 17290g 13 - 87672g 12 + 234696g n - 393054g 10 + 237158g 9 - 661318g 8 
+17401420/3g 7 + 338427338/3g 6 - 740356150g 5 - 2254648078/3g 4 + 47003116718/3g 3 
-125022788422/3g 2 + 34659398258g 

c [ 4) = -1632g 15 + 22400g 14 - 124488g 13 + 372512g 12 - 698512g n + 615952g 10 - 1348840g 9 
+3731152g 8 + 454575800/3<? 7 + 5824320064/15g 6 - 17623120168g 5 + 322919474240/3g 4 
-857689960864/3g 3 + 5197801864336/15g 2 -452601849976/3^ 

c [ 4) = -1938g 16 + 28560g 15 - 172690g 14 + 571402g 13 - 1191044g 12 + 1346952g n - 2483954g 10 
+2574458/3<? 9 + 163363530g 8 + 5936891456/3g 7 - 427652053304/15g 6 + 49753283730g 5 
+1888295056504/3g 4 - 11165804886700/3g 3 + 39268765896758/5g 2 - 5885730502990? 

c j 4) = -2280g 17 + 35904g 16 - 234720g 15 + 851200g 14 - 1957920? 13 + 2660560g 12 - 4434752g n 
-2180496g 10 + 458206840/3g 9 + 9171471680/3g 8 - 96727802776/15g 7 - 7627486265248/15g 6 
+15874774585472/3g 5 - 23613288071104g 4 + 809239034177296/15g 3 - 300645818028224/5g 2 
+24960425046264g 

C JJ = -2660g 18 + 44574g 17 - 313344g 16 + 1236200g 15 - 3116218g 14 + 14699672/3g 13 

- - 23502482/3g 12 - 13194884/3Q 11 + 369135052/3g 10 + 11095980830/3g 9 + 260264951278/9g 8 

- - 14453314408166/15g 7 + 218881556715946/45g 6 + 45588473570314/3g 5 - 2179954835300908/9g 4 
+14858869355291746/15g 3 - 81689671393093936/45g 2 + 3792619258364714/3g 

c j 4) = -3080g 19 + 54720g 18 - 411672g 17 + 1755776g 16 - 4819864g 15 + 25683904/3g 14 - 41145344/3g 13 
I Continued on next page| 
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TABLE IV: The coefficients ch> '(d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility X4( v > d) = 
En=o^(d)v n of the spin- 1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
the coordination number q = 2d. 

-12513184/3<? 12 + 229186600/3g n + 11910644656/3g 10 + 55759673624g 9 - 1666280612320/3^'' 

-199407718261024/15g 7 + 1180721878106528/5^ - 4913912212688992/3g 5 

+18295751469102064/3g 4 - 188910340455291256/15g 3 + 198430901982546976/15g 2 
-5313569529072792.7 

45 = -3542g 20 + 66500g 19 - 533178g 18 + 2445042? 17 - 7267702g 16 + 14344984g 15 
-23696940<? 14 + 3616100/3g 13 + 12996288g 12 + 58969981544/15g n + 75196713528g 10 
+232751024986g 9 - 1367841862391708/45g 8 + 1437869303984316/5g 7 
-10204868446463528/45<? 6 - 12244568432836332g 5 + 4116746384585661118/45g 4 
-1535126389650383104/5g 3 + 7658022498562556906/15g 2 - 336280716221910846g 



B. The HT expansion of x&(v;d) 



TABLE V: The coefficients c4 6) (d) of the HT series expansion 
in powers of v = tantrfC for the susceptibility x&{ v ]d) = 
X^ n =o oo c ™ 6 ^ (d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

c (6) - 16 
C W = i 36q 

ci 6) = 616g 2 - 376g 

cl 6) = 2016g 3 - 2912g 2 + 1016g 

c? 6) = 5376g 4 - 12768g 3 + 8592g 2 + 16q 

cp = 12432g 5 - 41664g 4 + 42584g 3 - 1568g 2 - 11784g 

cp = 25872g 6 - 112560g 5 + 157752g 4 - 33552g 3 - 48064g 2 - 16488g 

c w = 49632g 7 - 266112g 6 + 480816g 5 - 221888g 4 - 37848g 3 - 380032g 2 + 409096g 

ci 6) = 89232g 8 - 569184g 7 + 1270752g 6 - 944832g 5 + 292432g 4 - 1414576g 3 + 364752g 2 
+2177328g 

4 6) = 152152g 9 - 1125696g 8 + 3010224g 7 - 3132864g 6 + 1795080<? 5 - 7977904/3g 4 
-32649872/3g 3 + 112890736/3g 2 - 78703720/3g 

cS = 248248g 10 - 2090088g 9 + 6537168g 8 - 8795136g 7 + 6811560g 6 - 10886696/3g 5 

-139989920/3g 4 + 294356480/3g 3 + 283502264/3g 2 - 233678136.? 

c { S = 390208g n - 3683680<? 10 + 13226664g 9 - 21879744g 8 + 21066864g 7 - 5828736<? 6 

-106751176g 5 - 245869072/3g 4 + 1897397904g 3 - 11488470128/3g 2 + 2172800952g 

cjj = 594048g 12 - 6214208g n + 25234352g 10 - 49604016g 9 + 57229344g 8 - 17574592g 7 

-177848112g 6 - 2536571600/3g 5 + 5835827312g 4 - 8921212960/3g 3 - 23548883584g 2 

+30251052208? 

c{? = 879648g 13 - 10098816g 12 + 45813768g n - 104423264g 10 + 141116904g 9 - 64357856g 8 
| Continued on next page| 
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TABLE V: The coefficients <4 6) (d) of the HT series expansion 
in powers of v = tanhftT for the susceptibility X6( v 'id) = 
Sn=o oo Cn\d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

-228295552g 7 - 2298496240g 6 + 16344675328/3g 5 + 200049349136/3g 4 - 946502181448/3g 3 
+1442426211616/3g 2 - 236344803688(7 

of} = 1271328g 14 - 15890784g 13 + 79721096g 12 - 206842064g n + 321905760g 10 
-209254488g 9 - 189367680g 8 - 4356596032g 7 - 24934282288/3g 6 + 731828251760/3g 5 
-1884079942960/3g 4 - 2540774076440/3g 3 + 14787583499552/3g 2 - 4788568201288g 
eg = 1798464g 15 - 24310272g 14 + 133721952g 13 - 389311104g 12 + 688157624g n - 590354816g 10 
+91288088g 9 - 6810138432g 8 - 120995015840/3g 7 + 5693635677872/15g 6 + 4530618613376/3g 5 
-561 280804441 76/3<? 4 + 175625859309352/3g 3 - 1134460813568032/15g 2 + 101715603123320/3g 
eg = 2496144g 16 - 36279360g 15 + 217215936g 14 - 701494976g 13 + 1392033552g 12 
-1485565136g n + 949304208g 10 - 9440366096g 9 - 276020877632/3g 8 + 3692327961344/15g 7 
+40259238539312/5g 6 - 145655300075488/3g 5 + 117263299297904/3g 4 + 1915937192446272/5g 3 
-5668662187831232/5g 2 + 939412612077584g 

eg = 3405864g 17 - 52961664g 16 + 342995808g 15 - 1217241984g 14 + 2684716216g 13 
-3419593264g 12 + 3097994592g n - 12255126320g 10 - 161764858632g 9 - 1523657245328/5g 8 
+234835462900528/15g 7 + 7119519453424/15g 6 - 867980933559704g 5 + 76620462792781184/15g 4 
-193297122994892848/15g 3 + 74970873477425792/5g 2 - 6350764077866232g 

eg = 4576264g 18 - 75806808g 17 + 528159264g 16 - 2043636096g 15 + 4966825416g 14 
-22001068024/3g 13 + 23949277744/3g 12 - 47908848800/3g n - 735961413304/3g 10 
-6769987372664/5g 9 + 836907198003728/45g 8 + 1061178583645504/5g 7 
-123680489462001464/45g 6 + 128471706713417032/15g 5 + 557653654570467832/45g 4 
-2046057519029273912/15g 3 + 13722905419902781064/45g 2 - 683992228275224776/3g 
eg = 6063904g 19 - 106600032g 18 + 795192408g 17 - 3332563584g 16 + 8858492640g 15 
-44537998400/3g 14 + 55191538840/3g 13 - 68853616240/3g 12 - 1006105970464/3g n 
-2925741243056g 10 + 38101102223000/3g 9 + 2578706022827728/5g 8 - 34953325549594928/15g 7 
-469486142941733296/15g 6 + 1067819443550786792/3g 5 - 7740457369453181728/5g 4 
+17085072081639612176/5g 3 - 55533128504820921824/15g 2 + 1509351684072604536g 

45 = 7934080g 20 - 147517216g 19 + 1173245136g 18 - 5295282672g 17 + 15294767808g 16 
-28639450752g 15 + 39499504784g 14 - 116974101616/3g 13 - 420148419632g 12 
-74964569272096/15g n - 4879661509088g 10 + 3832420365724432/5g 9 
+33799538542615616/9g 8 - 649819887125975392/5g 7 + 36485867573612822272/45g 6 
-5101763740998382112/5g 5 - 89646139040158737760/9g 4 + 255768024175232201184/5g 3 
-1458409290472901242384/15g 2 + 67555026246347718672g 



The HT expansion of X8(v;d) 
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TABLE VI: The coefficients c^ 8) (d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility Xs( v ', d) = 
(s) 

^2 n= o oo Cn (d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

= -272 

cj 8) = -3968g 

4 s ) = -28160g 2 + 16064g 

cp = -135680g 3 + 18534 V - 63104g 

4 8) = -508640g 4 + 1152960g 3 - 774352g 2 + 39376g 

J 8) = -1595264g 5 + 5139200g 4 - 5262848g 3 + 626432g 2 + 1092480g 

cp = -4374656g 6 + 18399040g 5 - 25800320g 4 + 7123072g 3 + 7771104g 2 - 451712g 

4 8) = -10799360g 7 + 56229888g 6 - 101418880g 5 + 50302720g 4 + 20789760g 3 + 28718336g 2 
-48450944g 

4 8) = -24493040g 8 + 152280128g 7 - 338701088g 6 + 255076000g 5 - 7388208g 4 + 191935520g 3 
-212236432g 2 - 180336624g 

c ( 8) = -51800320g 9 + 374682880g 8 - 996706304g 7 + 1028397312g 6 - 349984384g 5 

+1895565824/3g 4 + 2459602432/3g 3 - 16331584256/3g 2 + 12754084736/3^ 

c j 8) = -103307776g 10 + 852508800g 9 - 2649851776g 8 + 3506874624g 7 - 2074473600g 6 

+4719123136/3g 5 + 24219243904/3g 4 - 83386258816/3g 3 + 12799466528/3g 2 + 28516712832^ 

= -195979264g n + 1817221120g 10 - 6480531200g 9 + 10519474944g 8 - 8601005312g 7 

+4050302720g 6 + 31146711936g 5 - 146794482688/3g 4 - 276543258112g 3 + 2167133336320/3g 2 
-442103339136g 

eg = -356075200g 12 + 3665247872g n - 14777517040g 10 + 28499934320g 9 - 29492323520g 8 
+39123812096/3g 7 + 248923396192/3g 6 + 38629503424g 5 - 4714014444080/3g 4 
+7509887769376/3g 3 + 6420043301440/3g 2 - 4809699889232g 

c [ 8) = -623054080g 13 + 7049482752g 12 - 31748068352g n + 71050636800g 10 - 88843430016g 9 
+141501441280/3g 8 + 520697373184/3g 7 + 1366919398144/3g 6 - 12055727532544/3g 5 
-6493984716800g 4 + 64312423603200g 3 - 335264513674496/3g 2 + 57772980342656g 

c [ 8) = -1054685440g 14 + 13009509760g 13 - 64792589056g 12 + 165243515008g n - 242500649248g 10 
+163985982336g 9 + 289429527488g 8 + 4584901860992/3g 7 - 16199812722112/3g 6 
-184263912218624/3g 5 + 822200412947456/3g 4 - 153864094174240g 3 - 2216382891425920/3g 2 
+922208891532608g 

c j 8) = -1733636608g 15 + 23152640000g 14 - 126442211584g 13 + 362278232576g 12 - 611273245440g n 
+1556604583424/3g 10 + 1003359499648/3g 9 + 3663324892416g 8 - 1635463424g 7 
-969658650984192/5g 6 + 705546067707904/3g 5 + 10957704133482496/3g 4 - 44381419834132480/3g 3 
+103569486797574912/5g 2 - 28828949928421504/3g 

eg = -2775830992g 16 + 39900474240g 15 - 237223499840g 14 + 754843908416g 13 - 1441747180848g 12 
+4473073125664/3g n - 69911950064/3g 10 + 22061233174384/3g 9 + 65837014833920/3g 8 
-384138748886400g 7 - 21910654023360544/15g 6 + 18929575971622144g 5 - 134319711356791312/3g 4 
-78918132289769312/3g 3 + 1092518218137459408/5g 2 - 210107883695368240g 

c [ 8) = -4340913280g 17 + 66818684416g 16 - 429792378880g 15 + 1504494886400g 14 - 3212833578240g 13 
+3936047152128g 12 - 1715109963776g n + 40032013286144/3g 10 + 222805972125440/3g 9 
-1579582584876032/3g 8 - 32879942504028672/5g 7 + 181099743084921344/5g 6 
+356398939022048512/3g 5 - 1318003726217876480g 4 + 19130142996335340032/5g 3 
-23490140191433877504/5g 2 + 2038219551714236544g 

c [ 8) = -6645196800g 18 + 109051063616g 17 - 754773335296g 16 + 2883607564800g 15 

-6815490996224g 14 + 28995994718080/3g 13 - 21216015570880/3g 12 + 69781088168704/3g 11 

| Continued on next page] 
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TABLE VI: The coefficients {d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility Xs( v ', d) = 
(s) 

^2 n= o oo Cn (d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

+174184899405984g 10 - 1227284831435200/3g 9 - 143067775267929728/9g 8 + 84830862839461504/15g 7 
+43761296741698393216/45g 6 - 5137578702684892480g 5 + 57115043408829771200/9g 4 
+314661275511421183136/15g 3 - 3143073756109544998336/45g 2 + 173286428915305918784/3^ 
eg = -9977510400g 19 + 173882649600g 18 - 1288851943808g 17 + 5338095665664g 16 
-13846023547392g 15 + 66947182003712/3g 14 - 64440145279744/3g 13 + 41372635858944g 12 
+340483735795968g n + 885961624852736/3g 10 - 84604254215282432/3g 9 - 14803804957801 1648g 8 
+37110053873309522176/15g 7 - 22502517113219316736/15g 6 - 253892478849624551680/3g 5 
+1436655318187681722368/3g 4 - 17346420241533257325056/15g 3 + 19533815567719217760256/15g 2 
-540509926226296804992g 

eg = -14718406560g 20 + 271460904000g 19 - 2145801037072g 18 + 9579095190288g 17 
-27071038685440g 16 + 48866734786048g 15 - 56625604724256g 14 + 236514693001600/3g 13 
+1767533919116336/3g 12 + 30408195156033568/15g n - 39513935844411504g 10 
-488739512196072000g 9 + 137400352859917843312/45g 8 + 112018250037801599872/3g 7 
-18735633799569055366448/45g 6 + 4151594372466586494272/3g 5 - 4171180399934389754912/45g 4 
-53306960347891458594336/5g 3 + 382608439113501754024976/15g 2 - 19084733342231658023472g 



D. The HT expansion of xio(w,d) 



TABLE VII: The coefficients c ( n '(d) of the HT series expan- 
sion in powers of v = t&nhK for the susceptibility Xw( v ', d) = 
^2 n= o oo Cn°\d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

cf> = 7936 

c y 0) = 176896g 

4 10) = 1805056g 2 - 990976g 

ci 10) = 11977856g 3 - 15817472g 2 + 5331456g 

cjj 10) = 59835776g 4 - 131649408g 3 + 88440832g 2 - 7113344g 

c| 10) = 243571328g 5 - 764301824g 4 + 784337664g 3 - 132473088g 2 - 130771200g 

e| 10) = 848335488g 6 - 3485722240g 5 + 4893727872g 4 - 1578948352g 3 - 1335948160g 2 + 306159872g 

c[ 10) = 2611834368g 7 - 13319674368g 6 + 24015744896g 5 - 12725984768g 4 - 5562137984g 3 
-1363634176g 2 + 7148102656g 

4 10) = 7273863168g 8 - 44390201856g 7 + 98556826368g 6 - 76197347072g 5 - 5531303424g 4 
-23187248384g 3 + 55450859776g 2 + 16088832768g 

4 10) = 18637718528g 9 - 132572504064g 8 + 351628333056g 7 - 365569735168g 6 + 83564687744g 5 
-362811788800/3g 4 - 13891630208/3g 3 + 2898002725888/3g 2 - 2566639813120/3g 

= 44506981376g 10 - 361758938112g 9 + 1120211169792g 8 - 1480274059264g 7 + 724243422848g 6 

-1363066769792/3g 5 - 4514271066752/3g 4 + 22869242557184/3g 3 

Continued on next page 
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TABLE VII: The coefficients <4 10) (d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility Xw( v > d) = 
Xm=o oo c « ^ {d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

-13265911387264/3g 2 - 4264749231360g 

c[\ 0) = 100054712576g n - 915089077760g 10 + 3249120416256g 9 - 5242545724416g 8 + 3817905779072g 7 
-1699281827840g 6 - 9196873191808g 5 + 81171677099264/3g 4 + 41908080196608g 3 
-498024823056896/3g 2 + 110610239423232g 

c{^ ] = 213450555136g 12 - 2169762639616g n + 8706410938368g 10 - 16654558212864g 9 
+15937418585088g 8 - 21327477729536/3g 7 - 104234607142528/3g 6 + 147562279746944/3g 5 
+1278777987933440/3g 4 - 1015433727966592g 3 + 18678581451392/3g 2 +969932534412928^ 

= 434934497536g 13 - 4864839662592g 12 + 21799836371456g n - 48343942806016g 10 
+57122586837888g 9 - 91561796736256/3g 8 - 298714168559488/3g 7 + 6166179994112/3g 6 
+1764241975937664g 5 - 2494450072249088/3g 4 - 46019039628909440/3g 3 
+94924943316071680/3g 2 - 17299176815162368g 

c{™ ] = 851007283456g 14 - 10386976407808g 13 + 51466387045632g 12 - 130025893949952g n 
+182747477058432g 10 - 122315166083968g 9 - 228199423744512g 8 - 954749235840512/3g 7 
+13606511889524864/3g 6 + 13727427846676864g 5 - 317929887704789632/3g 4 
+396942583026124160/3g 3 + 348318619589400320/3g 2 - 227789642969636608g 
c { ™ ] = 1606033313280g 15 - 21240894101504g 14 + 115402828579584g 13 - 327562821168128g 12 
+533995375623424g n - 443300013544960g 10 - 401964295244160g 9 - 1341116402467072g 8 
+23002048231205120/3g 7 + 1170205670425126912/15g 6 - 284186929185738880g 5 
-2073697908388501504/3g 4 + 13412366082401583104/3g 3 - 103692674008647256832/15g 2 
+9998443905803681024/3g 

4g 0) = 2934352776960g 16 - 41801425881600g 15 + 247247136328192g 14 - 779579603325952g 13 
+1447311237172736g 12 - 1455008875567104g n - 386203604227584g 10 - 11585115503070208/3g 9 
+5851443127592192g 8 + 3711330462883151104/15g 7 - 1646276274944524544/15g 6 
-7038770289879316480g 5 + 72320179932609938176/3g 4 - 193773243490858548224/15g 3 
-251833309601144910592/5g 2 + 60569364218487037184g 

cj 1 ^ = 5207010082560g 17 - 79484527964160g 16 + 508664327018496g 15 - 1764964370434048g 14 
+3679402474964736g 13 - 4372487916998656g 12 + 792697521304832g n - 9363099001265152g 10 
-43826578490797184/3g 9 + 8577550308943181056/15g 8 + 28219659363435557888/15g 7 
-375358120756778684416/15g 6 + 42682963471399987328/3g 5 + 1923822031694136063488/5g 4 
-20539071570523947907328/15g 3 + 8971213908296949751552/5g 2 - 799998504590970669312g 

c[f ] = 8998464034560g 18 - 146528416869120g 17 + 1009106692323840g 16 - 3822882507487232g 15 
+8848871709642496g 14 - 36520718842539776/3g 13 + 18743142681086720/3g 12 
-62988187629452800/3g n - 81387924415011456g 10 + 15571764869132276608/15g 9 
+383120271984801356288/45g 8 - 234112884730147594496/5g 7 - 13524450156382699305344/45g 6 
+12823863053571105585024/5g 5 - 249608017584109802434688/45g 4 

-23073655006552418847872/15g 3 + 893671518686969433181184/45g 2 - 567093417663341 19652096/3g 
c{l 0) = 15180194789760g 19 - 262646253492480g 18 + 1937297924305920g 17 - 7959422215602176g 16 
+20268674276510464g 15 - 95184859339027456/3g 14 + 74003803994501888/3g 13 
-139719129016333568/3g 12 - 242458510757969792g n + 4426762990994442496/3g 10 
+23718695776047253120g 9 - 375197918700176317696/15g 8 - 7237483465040656385536/5g 7 
+85047934197066602226944/15g 6 + 17815613791796507938176g 5 - 2607613209957224783747584/15g 4 
+7100994379431891533783168/15g 3 - 8373915013831095585411584/15g 2 
+236513875791114976989696g 

4J 0) = 25050011226240g 20 - 458902348544640g 19 + 3610245489776640g 18 - 15993366595831680g 17 

+44460160450820608g 16 - 78084468994640896g 15 + 78194195639666944g 14 

| Continued on next page] 
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TABLE VII: The coefficients c ( n '(d) of the HT series expan- 
sion in powers of v = tannic for the susceptibility Xw( v > d) = 
^2 n= o oo Cn°\d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

-320963319731515136/3<? 13 - 570503588763901312g 12 + 21341525437006926976/15<? n 
+51382154702027327744^ 10 + 757401336790171922432/5g 9 - 33169991145745310954752/9g : 
-11264809722306476971136/3g 7 + 8724305421769817690595712/45.? 6 
- 13644519855475631409950336/15g 5 + 1 1216450343505848520675328 /9q 4 
+11162831891410885080944768/5g 3 - 126554341938026732203131904/15g 2 
+6970666668081782114190720q 



E. The HT expansion of xi2(v;d) 



TABLE VIII: The coefficients & '(d) of the HT series 
expansion in powers of v = tanh.fr for the susceptibility 

(12) 

Xi2(v,d) = J2 n =o,oo °n (d)v n of the spin-1/2 Ising model on 
a (hyper)-simple-cubical lattice of dimensionality d as poly- 
nomials in q = 2d. 



= -353792 
4 12) = -11184128q 
4 12) = -154918400g 2 + 83051264g 
ct} 2) = -1351633920^ 3 + 1745895424g 2 - 585693184^ 
c? 2) = -8658773760<? 4 + 18669181440g 3 - 12553458432g 2 + 1241679616^ 

J 12) = -44288864256g 5 + 136435353600g 4 - 140287833088g 3 + 283016581 12g 2 + 19680018432g 
c j- 12) = -190551588864g 6 + 769929538560<? 5 - 1082384808960<? 4 + 386087227392g 3 
+264833080832g 2 - 88499933952g 

4 12) = -714269736960<? 7 + 3587327557632< ? 6 - 6469828485 120<r> + 3612002744320g 4 
+1473372426240g 3 - 380443199488g 2 - 1279808036864g 

4 12) = -2391664404480g 8 + 14393298296832g 7 - 31931620611072g 6 + 25342290355200g 5 

+2623725271552g 4 + 1090349782016g 3 - 14344234480128<? 2 - 840662029824g 

c ( 12) = -7287432355840<? 9 + 51178341826560<? 8 - 135515864580096g 7 + 142605665599488g 6 

-25221009321984g 5 + 49670513813504/3g 4 - 88730929958912/3.? 3 - 616452990285824/3g 2 

+626694374567936/3^ 

cj 1 ^ = -20494654415872g 10 + 164641056652800g 9 - 508612988841984g 8 + 674985995314176g 7 
-292637469406464g 6 + 319963052438272/3g 5 + 933634476506368/3g 4 - 6908597542308352/3g 3 
+5810927082614528/3g 2 + 712547619336960g 

c{\ 2) = -53792941017088<? n + 486701276938240g 10 - 1723075981086720<? 9 + 2779814633496576g 8 
-1875623051387904g 7 + 627651365216256g 6 + 3038947079074816g 5 - 36343321694500864/3g 4 
-4541280199228416g 3 + 136094192707704832/3g 2 -33234155640999936^ 

eg 2 ) = -132951019568640g 12 + 1338069948959744g n - 5351547829921280g 10 + 1020931229031 1680g 9 
-9262955488999680g 8 + 3645106505936896g 7 + 15429321176359936g 6 - 118905258348704768/3g 5 

-121040226141923840g 4 + 1218820247041000448/3g 3 - 150449878355645952g 2 

Continued on next page 
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TABLE VIII: The coefficients <4 12) '(d) of the HT series 
expansion in powers of v = tanh.fr for the susceptibility 
Xi2(v,d) = Sn=o,oo c ™ {d)v n of the spin-1/2 Ising model on 
a (hyper)-simple-cubical lattice of dimensionality d as poly- 
nomials in q = 2d. 

-234798668068850432Q 

= -311662549094400g 13 + 3454012141891584g 12 - 15422938056873984g n 
+34059565585510400g 10 - 38684966301057024g 9 + 19465049565298688g 8 + 56854453518541824g 7 
-85703458485491712g 6 - 2218501762470732800/3g 5 + 3619778929353328640/3g 4 
+12274347069953881088/3g 3 - 31957947017706407936/3g 2 + 6191270976593217536g 
c{ 12) = -697110544281600g 14 + 8436084083481600g 13 - 41644591004936192g 12 

+104685760328953856O 11 - 142679933981325312g 10 + 91914463877743104g 9 + 164968308753550848g 8 
-97434386494613504/ - 8225702434343863552/3g 6 - 5025700933230147328/3g 5 
+129303595089884297984/3g 4 - 222033723027268106240/3g 3 - 28411816242815580928/3g 2 
+68608468716171185408g 

= -1495254499614720g 15 + 19618972950528000g 14 - 106183823149867008g 13 
+299756354419372032g 12 - 476680111749890048g 11 + 383612063412850688g 10 
+373962329685188608g 9 + 133439709965193216g 8 - 21730914258414602240/3^ 
-431930711865523726336/15g 6 + 559493978162692917248/3g 5 + 163698341420861728768/3g 4 
-4707218563575137841152/3g 3 + 41124375194403225552896/15g 2 - 4136125039092031854592/3g 

c { ^ ] = -3088599584878080g 16 + 43673418631065600g 15 - 257320550871659520g 14 
+806768308656407552g 13 - 1467088196610935808g 12 + 1433944886486448128g n 
+556492490967948288g 10 + 1183366786353330176g 9 - 41636447435305676288/3g 8 
-415567921111926978560/3g 7 + 2011839579073394639872/5g 6 + 7952031547616848432640/3g 5 
-37897912866498410360320/3g 4 + 13077014203165422923776g 3 + 61638346384830787652608/5g 2 
-21409488905851651325440g 

c[\ 2) = -6166074284697600g 17 + 93474073271255040g 16 - 595872544132915200g 15 
+2055801777383424000g 14 - 4210415450621685760g 13 + 4873621078494318592g 12 
-220414900643438592g fl + 4510103641488453632g fo - 15880150474798704640g 9 
-450737986045982908416g 8 + 340344532355091195904/15g 7 + 221729576541278514829312/15g 6 
-32909314424528161202176g 5 - 351345139137285184575488/3g 4 + 8625784092226962208342016/15g 3 
-4070271793983798697129984/5g 2 + 374097575631435887130624g 

c{f ] = -11934548304806400g 18 + 193080915347523840g 17 - 1324578901007831040g 16 
+4989686958615091200g 15 - 11372156119571563008g 14 + 45791527096812011008/3g 13 
-17458751461905714688/3g 12 + 41937490887739925504/3g n + 31933338431806389760/3g 10 
-1153805055301468792576g 9 - 30373341839243953702400/9g 8 + 692581858601325041179904/15g 7 
+2590377436954334382597376/45g 6 - 39324075106561064944 70912/3g 5 
+33593400710271005522988032/9g 4 - 28199426216958322509515264/15g 3 
-293998710704524244730689536/45g 2 + 22806424112480971221711104/3g 
c[l 2) = -22455553019827200g 19 + 386163946273996800g 18 - 2837525699808184320g 17 
+11593375867440783360g 16 - 29118875151207942144g 15 + 133633187409078124544/3g 14 
-87651101589051000832/3g 13 + 122420035837644673024/3g 12 

+378885862990807923712/3g n - 2460433180576227344384g 10 - 469493202644201 14518016/3g 9 
+272796345569841313445888/3g 8 + 11043097785769283749154816/15g 7 
-76587116791527836322676736/15g 6 + 535138423569783282610176g 5 
+71151625301216987466575872g 4 - 1141571343228168064195402752/5g 3 
+4262414727355334828066523136/15g 2 - 123164411164904424887457792g 

c { ^ ] = -41170315007028480g 20 + 749910487155648000g 19 - 5877540461959246080g 18 
+25896572589293579520g 17 - 71102932389431354880g 16 + 122229074249888176128g 15 
-108969528004928226304g 14 + 357322559527178232832/3g 13 + 458200402892033175296g 12 
-66914070629597098638848/15g n - 47955624592773370635520g 10 + 80984110082112633531648g 9 

+134754839963129354103603712/45g 8 - 118043508845640078600391168/15g 7 
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TABLE VIII: The coefficients c ( n LZ) (d) of the HT series 
expansion in powers of v = tanh.fr for the susceptibility 
Xi2(v, d) = J2n=o,oo c ™ {d)v n of the spin- 1/2 Ising model on 
a (hyper)-simple-cubical lattice of dimensionality d as poly- 
nomials in q = 2d. 

-3966222753465191519663303168/45g 6 + 1733809313650138643157406208/3g 5 
-52123638124221675719579086592/45g 4 - 673442963130624721555058688/5g 3 
+49780845770293016030065316096/15g 2 - 3131394174087958521457421568g 



F. The HT expansion of xu(v,d) 



TABLE IX: The coefficients ch^ (d) of the HT series expan- 
sion in powers of v = tanh.fr for the susceptibility Xu( v i d) = 
Sn=o oo Cn A \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

ci 14) = 22368256 
C J 14) = 951878656<? 

4 14) = 17171485696g 2 - 9061682176g 

cl 14) = 190346960896g 3 - 242153541632g 2 + 81035140096g 

c| 14) = 1519012888576g 4 - 3228372738048g 3 + 2173094684672g 2 - 241910319104g 

4 14) = 9522179620864g 5 - 28942242217984g 4 + 29810692870144g 3 - 6697031569408g 2 
-3638385785856g 

c { 6 u) = 49521284890624g 6 - 197607542865920g 5 + 278159316692992g 4 - 106288041009152g 3 
-61351205630976<? 2 + 25207593720832^ 

4 14) = 221721122013184g 7 - 1100726949249024g 6 + 1986223294693376g 5 - 1154500961435648<? 4 
-426102323097600g 3 + 245451202273280g 2 + 271935058287616^ 

4 14) = 877590034038784g 8 - 5224913761435648g 7 + 11588424057585664g 6 - 94104191273861 12q 5 

-1000744824053760g 4 + 1075295221833728<? 3 + 4029547520798720<? 2 - 338914506043392(7 

4 14) = 3132002019784704g 9 - 217768714517381 12g 8 + 57607925238919168<? 7 - 61388804936712192g 6 

+9491513408580608g 5 + 6447340566738944/3g 4 + 44754385667900416/3g 3 + 154173876921395200/3g 2 

-182161361411537920/3g 

4o 4) = 10232670234040320g 10 - 81443029784979456g 9 + 251210594330341376<? 8 

-335497615196946432g 7 + 135461448552674304g 6 - 35792119632002048/3g 5 - 222702217883454464/3g 4 
+2362461652009822208/3g 3 - 2456278422725023744/3<? 2 - 111726801462610944^ 

c[\ A) = 30973029286748160g n - 277822048171008000<? 10 + 981629865368946688g 9 

-1587563073560690688g 8 + 1020762474263646208g 7 - 209702033422557184g 6 - 1173173825248046080g 5 
+16683124950400256000/3g 4 - 1176998384765688832g 3 - 43334277608785086464/3g 2 
+11767439556985869312g 

eg 4 ) = 87695613221928960g 12 - 875516256400220160g n + 3493448597460332544<? 10 
-6664923981790646272<? 9 + 5837821311797616640g 8 - 5668522943666667520/3g 7 
-22994941740443248640/3g 6 + 25365453641949280256g 5 + 103451675164792563712/3g 4 

-526401211985641717760/3g 3 + 312941970439251582976/3g 2 + 65165416306793678848g 
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TABLE IX: The coefficients On A) (d) of the HT series expan- 
sion in powers of v = tanhfT for the susceptibility Xuiv, d) = 
X^ n =o oo Cn A \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

c£ 4) = 234088827521064960g 13 - 2574836787374161920g 12 + 11467569791287357440g n 
-25287282865045250048g 10 + 279561654850245 73440g 9 - 38557986919775731712/3g 8 
-104508826006886752256/3g 7 + 255375882161599811584/3g 6 + 952361429739659500544/3g 5 
-846995991975855159296g 4 - 1139738243810908156928g 3 + 12549696909467484780544/3g 2 
-2597769507194729141248g 

c^ 4 4) = 592949907428597760g 14 - 7125199541722030080g 13 + 35076261083993886720g 12 
-87958340403507830784g n + 117329974060356933632g 10 - 71843845823446751232g 9 
-122586910874966786048g 8 + 669004110698680303616/3g 7 + 4744396665028853771264/3g 6 
-3669495689471793658880/3g 5 - 57018764157647422594048/3g 4 + 41105792891731103113216g 3 
-25186090822497939331072/3g 2 - 24319035703272885756928g 

c[ 1A) = 1433048994515435520g 15 - 18678993309190717440g 14 + 100804343240191549440g 13 
-283697078951770030080g 12 + 443223271221351739392g n - 1033637533951132155904/3g 10 
-1011217207474809448448/3g 9 + 460483675184064487424g 8 + 5565906744552960131072g 7 
+43919638020575731011584/5g 6 - 347024996061061408676864/3g 5 + 239690733943793419577344/3g 4 
+1866879933655208536607744/3g 3 - 6335397272764050983266304/5g 2 
+2000581556284493794393088 /3q 

c[ lA) = 3319867725632901120g 16 - 46653317999297740800g 15 + 274061753900784107520g 14 
-856289832024370053120g 13 + 1533976830053824094208g 12 - 4377264783509685641216/3g n 
-1930294775907870310400/3g 10 + 2043781452722642796544/3g 9 + 44749069561807159070720/3g 8 
+368688640717095980146688/5g 7 - 6020966558926158905196544/15g 6 - 959320791731322801360896g 5 
+21141616501494931472162816/3g 4 - 144951016231397003748696064/15g 3 
-11761559674844284219748352/5g 2 + 8971329699132498653917184g 

c[ lA) = 7401483245876689920g 17 - 111549432110458552320g 16 + 708955002035236392960g 15 
-2436938752035708026880g 14 + 4926616716589229905920g 13 - 5571945225556920684544g 12 
-164482664925404280832g fl + 542451596580084125696/3g 10 + 90500305739342358821888/3g 9 
+4935275039197882432403456/15g 8 - 3828311185429464106921984/5g 7 
-43083052926904186971451392/5g 6 + 91947713254873005894375424/3g 5 
+150183300690139346326908928/5g 4 - 1386080926679917894740688896/5g 3 
+2150316704653739855171315712/5g 2 - 204220289203880776017288192g 

c[ 1A) = 15934593393411025920g 18 - 256381800437179069440g 17 + 1753506681434233835520g 16 
-6580315733798396805120g 15 + 14826378622802061465600g 14 - 19494270830426973841408g 13 
+5840961809226396641280g 12 - 12379160817027722752000/3g n + 114757699742723572333568/3g 10 
+16488703030550456426288128/15g 9 + 9566178324343212097599488/45g 8 
-191473112196390312892758016/5g 7 + 1539992983726049687744905216/45g 6 
+10601363112007581786085537792/15g 5 - 114323904676714821742061170688/45g 4 
+34146144587643522501476139008/15g 3 + 103223347344087124949760379904/45g 2 
-10901839332941488863023727616/3g 

c{\ A) = 33226125696362987520g 19 - 568427635729473269760g 18 + 4164157054396259312640g 17 
-16947979998253087703040g 16 + 42132111162172651745280g 15 - 63257522133446069288960g 14 
+36355926368650676676608g 13 - 74447257194277928488960/3g 12 - 54484321215182590646272/3g 11 
+9097409010673707965517824/3g 10 + 24262474856816972871058432/3g 9 
-1714006446732822316017821696/15g 8 - 4665520708535731448468758528/15g 7 
+61022631460493300164426153984/15g 6 - 5972682674402348291090053120g 5 
-468681727160463278326158270464/15g 4 + 1905595883401702552479227158528/15g 3 
-2524382647775783152293231110144/15g 2 + 7477639870846364483 1068077056g 

4n 4) = 67276931969120501760g 20 - 1219446678396708372480g 19 + 9528887956728128378880g 18 
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TABLE IX: The coefficients c { n ' (d) of the HT series expan- 
sion in powers of v = tanhfT for the susceptibility Xuiv, d) = 
X^ n =o oo Cn A \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

-41822293160168382044160g 17 + 113760507839982661263360g 16 - 192236223612750884290560g 15 
+156489073462191969669120g 14 - 110241775369845605625856g 13 -919955615111802494734336/3? 
+36119687323866478101979136/5g n + 38236807443309719843198976g 10 
-1191331056719902439094376448/5<? 9 - 19547313607379395121578246144/9g 8 
+180510166791180821800942637056/15g 7 + 1605759665205848894990726950912/45g 6 
-5794973282407102465659826425856/15<? 5 + 8663267176107152845551619469312/9g 4 
-2296147211477278179472071303168/5g 3 - 22163092582660818160691985946624/15g 2 
+1673593062235423907376589068288Q 



G. The HT expansion of Xie{v;d) 



TABLE X: The coefficients d 16) (d) of the HT series expan- 
sion in powers of v = tanhi^ for the susceptibility Xi6( v > d) = 
X^ n= o oo Cn 6 \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

A 16) = -1903757312 

c | 16) = -104932671488g 

ci 16) = -2389096202240g 2 + 1247014436864^ 

c| 16) = -32769353973760g 3 + 41234501730304g 2 - 1377981 1467264g 

4 16) = -318434742599680g 4 + 669638839173120g 3 - 451194468442112g 2 + 54187912339456g 

4 16) = -2398328708005888g 5 + 7216293157273600g 4 - 7443684913840128g 3 
+1796451694804992g 2 + 810471469645824g 

4 16) = -14814711540744192g 6 +58553079162798080g 5 - 82516960690176000g 4 

+33131931010138112g 3 + 16568564934164480<? 2 - 7826521298747392g 

c (i6) = _ 780 oil93991536640g 7 + 383758520681496576g 6 - 692893232058204160g 5 

+415590082033745920g 4 + 138050755428352000g 3 - 112563191669850112g 2 
-67250352764125184^ 

c ( 16) = -359881375772160000g 8 + 2124537375912247296g 7 - 4712048562925756416g 6 
+3901175774945075200g 5 + 382668819973588992g 4 - 830697243533570048g 3 
-1253425152787441664g 2 +246687187735013376^ 

c ( 16) = -1485475101516595200g 9 + 10246573574686310400g 8 - 27091221120030670848g 7 
+29207910483003244544g 6 - 4295104640277151744g 5 - 12340496168646148096/3<? 4 
-19033354976936984576/3g 3 - 44592692556107743232/3g 2 + 62059824450177499136/3^ 
c[ l 6) = -5573914967367352320g 10 + 44032574907530035200g 9 - 135681242051050831872g 8 
+182431884556188581888g 7 - 70873489496818745344<? 6 - 32360324729745539072/3g 5 
+67706322643199107072/3g 4 - 919115973967527559168/3g 3 + 1110469663689548767232/3g 2 
+4028078481045454848g 

c\\^ = -19254430629988270080g 11 + 171497359126619750400g 10 - 605111828654049361920g 9 

| Continued on next page | 
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TABLE X: The coefficients c4 16) (d) of the HT series expan- 
sion in powers of v = tanhiT for the susceptibility Xi6( v ', d) = 
X^ n =o oo Cn 6 \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

+982109969155263234048g 8 - 613201437461336227840g 7 + 52147748159124078592g 6 
+540787451457111064576g 5 - 8293323481894075432960/3g 4 + 1601119273086217224192g 3 
+15799373824615516340224/3g 2 - 4836104293955961716736? 

cj 16) = -61858597238822092800g 12 + 613495164956018442240g n - 2443811528379954585600g 10 
+4669086468011990999040g 9 - 3997355520769368145920g 8 + 3048460248858627997696/3g 7 
+13023044456997156921344/3g 6 - 48507816777292733022208/3g 5 - 24805283810470472028160/3g 4 
+83595989933429336940544g 3 - 191790280034196466143232/3g 2 - 19678447759486369759232.? 

c[f ] = -186378077798827622400g 13 + 2037305764689922621440g 12 - 9056133528765221437440g u 
+19970129771694656716800g 10 - 21679251241263612198912g 9 + 26990843334412000231424/3g 8 
+70554377277375853297664/3g 7 - 213211678483330007105536/3g 6 - 139778987713044081672192g 5 
+1664707021142943756255232/3g 4 + 839278592660028565356544/3g 3 
-5660286687163968965574656/3g 2 + 1258500333536581549654016^ 

c [ 16 ) = -530309836482532147200g 14 + 6335138766411898060800g 13 - 31121484056466262917120g 12 
+77967410300989969858560g n - 102461246405166576254976g 10 + 59417543428498750226432g 9 
+97390176985860741332992g 8 - 762964738187005595942912/3g 7 - 2771862731613706912448512/3g 6 
+1762940830895704098471936g 5 + 27268466707857705964199936/3g 4 
-73144748403820602006495232/3g 3 + 29494003189156796786991104/3g 2 
+9859878527425768111505408^ 

c [ 16) = -1433313857854681251840g 15 + 18579093268881997824000g 14 - 100041765469328009134080g 13 
+281099599114713305579520g 12 - 433684248549726146396160c/ 11 + 325666378885024346406912g 10 
+314189928509364533886976g 9 - 778640166888945468899328c; 8 - 12345110484432923427930112/3g 7 
-6029605745496073616162816/15g 6 + 74183532794057737351725056g 5 
-299336476697102553835700224/3g 4 - 812355558976921220939382784/3g 3 
+10079179148554636101429231616/15g 2 - 1113855886758632875032739840/3g 

c [ 16) = -3698167731092948367360g 16 + 51697759510997113651200g 15 - 302989963080538533273600g 14 
+944731356246521622405120g 13 - 1674213903991364760268800c; 12 + 1552712554992677806497792c/ 1 
+720447356096041891672064g 10 - 6282927250512377689477120/3g 9 - 13916774158254107683837952g 8 
-109261731058784289328144384/3g 7 + 5247550300007197475660517376/15g 6 
+255075570728706442667266048g 5 - 12816733701499485987675908096/3g 4 
+21265133695803460802539581440/3g 3 - 2749516121165670024115259392/5g 2 
-4346747332669088027247523840g 

C J 16) = -9147798671639519232000g 17 + 137186452951627556782080g 16 - 869813208451446113894400g 15 
+2982818202880574005248000g 14 - 5973178395596154214809600g 13 + 6619544164878738447728640g 12 
+481000794676112127098880g 11 - 4899603407702683328774144g 10 - 110565361480733286376570880/3g 9 
-704983319085229368702533632/3g 8 + 16103452583665361645089193984/15g 7 
+75829256330761417932755566592/15g 6 - 77977148669047737921882554368/3g 5 
+1811968389847140769282719744g 4 + 2257670174798245419856825286656/15g 3 
-1302173999037493066633818865664/5g 2 + 128088962284953555125162901504g 

c ( 16 ) = -21773726436257454489600g 18 + 348689272216526236385280c/ 7 - 2379062154498192511303680g 16 
+8904943167941543976960000g 15 - 19894429006009620467220480g 14 + 25707317955109416860549120g 13 
-6350952062228959616368640g 12 - 8916914984431825074126848g fl - 219363621455006506598907904/3g 10 
-3025488708891444056995053568 /3g 9 + 16821728067680996475046395904 /9g 8 
+464207279270231673467090960384/15g 7 - 3281882880491767934129981210624/45g 6 
-1193128885900196613472098156544/3g 5 + 16617460149856001792874156007424/9g 4 
-32155737405107925626486996885504/15g 3 - 32556846223405514512001694785536/45g 2 

+6046576182698989570541898727424/3g 
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TABLE X: The coefficients d '(d) of the HT series expan- 
sion in powers of v = tanhK for the susceptibility Xw(v; d) = 
X^n=o oo Cn 6 \d)v n of the spin-1/2 Ising model on a (hyper)- 
simple-cubical lattice of dimensionality d as polynomials in 
q = 2d. 

c Jg 6) = -50030037863288065228800g 19 + 852101290521418412851200<? 18 - 6227047101542461421322240g 
+25275614189702156621905920g 16 - 62352271739853593258557440a 15 + 92187604552840911877242880^ 
-47835588891081063187742720g 13 - 45489007942380861 13460224g f2 
-237947979009622578034180096/3Q 11 - 10366035347757135714664644608/3g 10 
-1265874199130683595892228096g 9 + 123225127752271144007575142400g 8 
+634660532275614856202776477696/15g 7 - 48783751300574925991281498914816/15<? 6 
+23362355833110189429704435826688/3<j 5 + 40370908413300708220306030919680/3<j 4 
-1202743242344609002546580936851456/15g 3 + 1712098811966697109750802892783616/15g 2 
-52084653311579531984812997640192^ 

cgf) = -111284638050657492664320g 20 + 2008604722638854258688000g 19 
-15657114688558678873128960*? 18 + 68528531116644896241008640g 17 
-185091246922911993627525120g 16 + 308456737644166806503424000a 15 
-233888851657298026447503360g 14 + 67391655489588897810415616a 13 
+141274859918430596363763712g 12 - 50591878074073047993295863808/5O 11 
-24823899558370202782000799744g 10 + 362123860232776118591214739456g 9 
+64086239757063227627779602399232 /45g 8 - 207863878816724357223254796271616/15g 7 
-283156467508485410552545878499328/45g 6 + 828660290657745177421922636824576/3g 5 
-37071865469327969561767228813586432/45g 4 + 660261762816077560188055602307072g 3 
+10536441619319090528029264686166016/15g 2 - 1039099543218814255033842395971584g 
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